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Cosmological models involving an interaction between dark matter and dark energy have been 
proposed in order to solve the so-called coincidence problem. Different forms of coupling have been 
studied, but there have been claims that observational data seem to narrow (some of) them down 
to something annoyingly close to the ACDM model, thus greatly reducing their ability to deal with 
the problem in the first place. The smallness problem of the initial energy density of dark energy 
has also been a target of cosmological models in recent years. Making use of a moderately general 
coupling scheme, this paper aims to unite these different approaches and shed some light as to 
whether this class of models has any true perspective in suppressing the aforementioned issues that 
plague our current understanding of the universe, in a quantitative and unambiguous way. 



I. INTRODUCTION 

Several sets of observational data such as Type la su- 
pernovae (SNIa), Cosmic Microwave Background (CMB) 
and large scale structure surveys, when combined, in- 
dicate that we live in a nearly flat, low-matter-density 
(Omo ^ 0-3) universe whose expansion is speeding up 
at present. The driving source of this cosmic accelera- 
tion is nevertheless still poorly understood. The simplest 
candidate is Einstein's cosmological constant A, assumed 
to stand for the vacuum energy density (pA — A/SttG). 
However its tiny value, as inferred from observations, is 
some 50 (120) orders of magnitude below conservative 
(aggressive) estimates given by quantum field theory and 
explaining this discrepancy is one of the most difficult 
problems in theoretical physics. This constitutes what 
is sometimes called the cosmological constant problem 
(CCP), whose solution probably require the discovery of 
a yet unknown underlying symmetry capable of inducing 
a nearly perfect cancellation of the vacuum energy den- 
sity value (candidates such as supersymmetry do not help 
since they are broken at high energy scales) or at least of 
its contribution to gravity. Such fine-tuned cancellation 
is thought to be very unlikely, and the problem is usually 
addressed by what is considered a more reasonable hy- 
pothesis: a complete suppression of the vacuum energy 
density contribution to the gravitational sector. This in 
turn requires a different explanation for the aforemen- 
tioned cosmological observations. 

Possible candidates discussed in the literature usually 
fall into two categories: exotic components with negative 
pressure (dark energy or simply DE) or proper modi- 
fications of general relativity which become relevant at 
cosmological scales. In either case different observational 
and theoretical problems arise, and this made the soil 
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fertile for a plethora of models to bloom. On the the- 
oretical side we identify three possible issues which we 
will dub: the coincidence problem (CP), the DE energy 
density initial condition problem (DEICP) and the sen- 
sitivity to initial conditions (SIC) . These can be strongly 
related (as in the case of any A-like dark energy model, 
defined below) or even dependent on the physical inter- 
pretation underlying a particular model (as opposed to 
being a problem of the model itself), as we will clarify 
below. Another important point is that all three prob- 
lems often impose fine-tuning to observationally viable 
models. 

The CP can be stated as: why only re- 
cently (in terms of redshift) has the DE en- 
ergy density become comparable to the matter one, 
since both components are usually assumed inde- 
pendent and thus scale in different ways? The 
DEICP appears when one interprets dark energy as a 
proper field, but not when one works in modified gravity 
theories. In other words: why is the initial (by which 
we mean just after inflation, z ~ 10^^) value of the DE 
energy density much smaller than what we would expect 
from equipartition considerations? It seems natural to 
expect that after inflation the different fields in nature 
would have energy densities with the same order of mag- 
nitude. Thus a reasonable value for the ratio of energy 
densities between DE and radiation, pbe/Pi- , would be 
~ 10~^ or 10"'^. The SIC is a measure of the robustness 
of a model to different initial conditions, which usually 
can be thought in terms of basins of attraction. Models 
with a larger SIC require a more fine-tuned initial condi- 
tion than those with a smaller one. The amount of SIC is 
not always considered an issue, as the initial condition for 
radiation itself can be thought of as fine-tuned in order to 
match present measurements of the CMB temperature. 
Nevertheless, models which allow a broader range of ini- 
tial conditions are usually preferred over those which put 
stringent constraints on such a range. The relevance of 
each of these three issues (CP, DEICP and SIC) can be 
disputed, but our main intention here is to clarify their 
distinction, as is not uncommon for ambiguity among 
them (and even with the CCP) to arise in the literature. 
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Before we proceed, a remark on notation is in order. 
To make explicit our phenomenological approach (which 
includes for instance modified gravity models) in our dis- 
cussions we will denote by the indexes "x' and "de" all 
physical quantities related to DE, but we will reserve the 
latter exclusively to the cases in which only the field in- 
terpretation makes sense. Furthermore, we will use the 
index suffixes "o" and "i" to denote respectively quan- 
tities evaluated today and at the end of inflation. On 
an altogether different issue, we will refer to la, 2a and 
3a as a shorthand notation for 68.3%, 95.4% and 99.73% 
confidence levels, respectively, even though this would 
only be rigorously the case for gaussian likelihoods. 

In ACDM (henceforth ACDM will denote any model 
which behaves like a cosmological constant but which is 
not described by the vacuum energy density) or in any 
other constant-uij; CDM model, the CP and DEICP can- 
not be jointly addressed, and either one requires a fine- 
tuning of a parameter (e.g. pDEi)- To have a A domi- 
nated universe today, the initial DE energy density has 
to be set dozens of orders of magnitudes smaller than 
the initial matter energy density. However this may not 
be the case, for instance, for quintessence tracker scalar 
field models [T]. Attractor-like solutions appear in this 
case such that, for a wide range of initial conditions, 
the field eventually reaches a definite cosmic evolutionary 
track at which f2mO ~ ^^deo- Assuming that the back- 
ground energy density decreases with the scale factor as 
Pback oc a"", we call tracker those solutions for which the 
energy density of the tracker field scales as px oc a~'", 
where m depends on the equation of state of the back- 
ground but is such that m < n. Scaling, on the other 
hand, occurs when the ratio Px/pm is constant [5], that 
is, when m ^ n. Adjusting the onset of tracking is one of 
the difficulties of quintessence tracker solutions. For in- 
stance, starting from equipartition conditions after infla- 
tion it was shown in [T] that, in inverse power law tracking 
quintessence models {V{(j>) cx </>""), only for a > 5 will 
the fleld start tracking before matter-radiation equality. 
However current observations impose a < 1 [21 H], and 
in this case the fleld would only start tracking at recent 
epochs. Since DE must develop negative pressure around 
z 1 in order to drive acceleration, this demands a flne- 
tuning to the fleld potential. Therefore, although the 
DEICP can be addressed by these models, they pose no 
clear advantage as far as the CP is concerned, as the lat- 
ter could be reformulated as: why do we live in a special 
epoch when the fleld just started tracking? In order to 
circumvent this difficulty, models with a non-canonical 
kinetic term (k-essence) O [6] have been proposed as an 
alternative. In k-essence, the fleld reaches a scaling so- 
lution during the radiation dominated era, but dynamics 
(rather than an adjustment of a parameter) triggers a 
transition when non-relativistic matter starts to domi- 
nate, after which k-essence follows an accelerated attrac- 
tor. However, a good candidate with a large basin of 
attraction is still lacking [7]. 

Here we shall follow the approach of [8 , and address 



the CP by requiring a duo-scaling cosmology. We seek 
models whose phase space exhibit (at least) two flxed 
points: one responsible for the matter dominated era 
(henceforth MDE) and the other for the present dark 
energy dominated acceleration. The former needs to be 
a saddle point and the latter (preferably) an attractor. 
The CP could be solved by a model in which our present 
universe has reached (by which we mean it is close enough 
for a given criterium) such an accelerated attractor, as 
this would mean that the current cosmological energy 
distribution is not a transient phase but rather an un- 
avoidable and permanent regime. The existence of such 
stable and accelerated flxed point requires a coupling be- 
tween matter and dark energy |^, as otherwise the DE 
energy density would have to decrease with a~'^. 



In order to quantify both the DEICP and SIC, 
we will deflne two quantities: C = Pu^/Pri and 
A EE (p"f - p™'")/ Ip'xt"^] , where pr stands for the radi- 
ation energy density and P™/^^'" denotes the maximum 
and minimum values of Pxi/BEi that evolve to the present 
day observed values within la. Models which have small 
SIC are characterized by large values of A and vice-versa. 
Note that our choice of variable for measuring the SIC is 
a ratio of the size of the permissible region to the maxi- 
mum value of that range, and thus is not directly related 
with the proper size of the allowed region of initial con- 
ditions. In other words, A is a relative measure, not an 
absolute one. One is sometimes more interested in eval- 
uating the absolute width of the initial condition range 
and in this case {PxT^ — P'xi'^) could be used instead. One 
should bear in mind that newer and improved observa- 
tions will probably narrow our present uncertainties on 
Px and this could reflect on the initial range A, which 
should depend only on the model and not on the quality 
of our observations. This is avoided by computing the 
ratio A/Aacdm- Furthermore, our deflnition of A is not 
a good one in the cases where p™'*'^ is very close to zero 
and p""" is not. Finally, for the CP we will use the ratio 
Rz between DE and dark matter energy densities, where 
z stands for the redshift, as a measure of how close we are 
to the accelerated attractor. More speciflcally, a solution 
to the CP requires ~ Rq. 

Making use of a three-parameter model with a cou- 
pling scheme that generalizes some previous ones in the 
literature (section |ll]) , we analyse the phase space (sec- 
tion HI I and apply two different cosmological tests (sec- 
tion ^Vf: SNIa, as given by a combined catalog [TO] . 
and the so-called CMB shift parameter, as inferred from 
WMAP3 [11 . We look at three aspects: what is the 
parametric region that allows such duo-scaling cosmol- 
ogy (section IV ) ; how does the model cope with the DE- 



ICP inside that parametric region; and what is the SIC 
(section VI L We also propose a two-parameter toy model 



VII| for which there exists analytic solutions and 

in some cases. 



(section 

which has small SIC, alleviate the CP and 
also the DEICP. Finally, we note that by choosing freely 
all three parameters we can actually solve the CP. 
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II. THE DARK INTERACTIONS MODEL 

We consider the universe as filled by four components: 
radiation, baryons, and two coupled barotropic (latu 
sensu) fluids. One, a cold, pressureless, dark matter 
(CDM) and the other a negative pressure dark energy. 
We will denote them, respectively, by the subindexes r, 
b, c and x. For the sake of simplicity, our analysis will 
be restricted to models with constant equation of state 
parameter = Px/Px- Local gravity constraints limit 
greatly any possible interactions between dark energy 
and baryons. To allow such a coupling requires a mecha- 
nism to either make these constraints only effective in the 
present and not in the past, or to make the interaction 
range short enough, as achieved in the chameleon scalar 
field models [12l[13]. We therefore neglect this possibihty 
and focus on a coupling with dark matter alone. An in- 
teraction with radiation is also discarded on the grounds 
that any such reaction would not affect the dynamics of 
the system near the sought two scaling regimes, required 
to address the CP. On the other hand, a coupling with 
radiation might be desirable if one wants to address also 
the DEICP. Since this would introduce many difficulties 
of its own, we will not consider it further in this work. 

Our interest is in cosmological scaling solutions in 
a spatially flat Friedmann-Lemaitre-Robertson- Walker 
(FLRW) background metric with a scale factor a{t): 



ds^ = -d<2 + a2(t)dx2 



(1) 



The Friedmann equation in General Relativity is given 

by 



3i?2 = Mp^{p^ +Pc + Pb + Pr) 



(2) 



where Mp — SnG with G being the gravitational con- 
stant, and where the different pi make up the energy den- 
sity of the universe. In what follows we shall set Mp — 1. 
Different forms of coupling have been considered in the 
literature. We may take the coupling between CDM and 
DE to be such that 



dpx 
dN 



+ 3(1 -I- Wj;)px = -3Q{pc, Px), 



dpc 
dN 



(3) 



+ 3pc = iQiPcPx), 



where = In a, sometimes called the number of e-folds. 
Then, some of the (non- mutually exclusive) proposed 
forms for Q{pc,Px) so far are: 



QiPcPx 



(A - w^) 
A ^' 



PxPc 



Px 



VPx + Pc 
A {px + Pc) 
XPc 

, '^l Px Pc 



mils] 

m 

m 
[niiis] 

m 



(4) 



where, in addition to Wx, we take A and the different 
Xi to be constants. Note that positive values of Q{pc, Px) 
indicate that energy is being transferred from dark energy 
to dark matter, meaning that the latter will dilute slower 
than in the case without interactions. 

Our analysis will be restricted to the case"'^ 



Q{Pc, Px) = X^px + XcPc 



(5) 



which is a more general form than those found in |17j 
or 19J, but reduces to those cases, respectively, when- 
ever Xx = Xc = X and A^; = 0, Ac = A. It is, however, in 
a different class than the ones used in many interacting 
scalar field models. In particular, it is not equivalent to 
the coupling used in ^ and therefore it might circumvent 
the tough challenges regarding the solution of the CP im- 
posed there. Furthermore, in what follows we will always 
consider Wx to be a negative constant, whose value is the 
third free parameter of our model. 

One of the reasons for considering the above coupling 
scheme is that it is a straightforward extension of some of 
the other forms encountered in the literature. Nonethe- 
less, it has some distinctive features that the reader 
should bear in mind. First, when one introduces a cou- 
pling term which is proportional to Pc but not, say, PcPxi 
one has to be careful about what one means by "dark 
matter" . The reason is that even in the absence of dark 
energy, the former has a non zero equation of state pa- 
rameter. In fact, for p^. = we have — —Ac, as 
can be easily seen from eq. [3] Second, such coupling nat- 
urally allows for px and pc to become negative as the 
universe evolves, and the way to deal with this is dis- 
cussed in the following section. Third, the overall sign 
of Q{pc, Px) may change as the energy densities dilute. 
A late-time scaling leading to cosmological acceleration 
enforces a positive value today, but leaves open the pos- 
sibility of an opposite energy transfer in the past, which 
in turn could lead to a stronger formation of structures 
during the MDE. 

Rewriting eqs. (|3l using eq. (Isl we get 



dpx 
dN 
dpc 
dN 



+ 3(1 +Wx + Xx)px = 

-1-3(1 - Xc)pc = 3XxPx 



-3 Ac Pc, 



(6) 



This system does not admit analytic solutions in the 
general case^ but there are noteworthy exceptions for 
some particular ones, namely: (i) A^; = 0; (ii) Ac = 0; 
(iii) Ac = 1 and (iv) Wx — —l — Xx- For each of these 
possibilities, it is possible to decouple the above equa- 
tions. Case (i) is the one studied in fT^. The second 
case, while in principle worth investigating, turns up to 



^ It has only recently come to our attention that this proposed 

form was previously considered by 1211 . 
^ General solutions exist if one neglects baryons and radiation, as 

has been shown in |22| . 
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be incapable of alleviating the initial condition problem, 



as we will show in section VI Case (iii) is completely 



ruled out by observations, as will become clear in sec- 
tion |Vj The last possibility, however, is an interesting 
one. It allows simultaneously Wx and A^; to be close to 
— 1 and 0, respectively, while still exhibiting a distinct 
behaviour compared to ACDM. Not only can it easily 
satisfy the observational tests we employed, but it also 
addresses both the CP and DEICP. We will return to this 
particular choice of parameters in section |VII[ 



III. FIXED POINTS 

In order to study the dynamics of the model, we shall 
introduce the following set of variables, 

X = J-^ Y = — ^ 
i/2 3 ' il2 3 ' 



z = 




1 



Pr 
3 ' 



(7) 



and rewrite the continuity equations and eq. Q as 
X [3wx{X -l) + u^- 3Aj - 3Ac Y, 



= Y [3wxX + u^ + 3Xc] + 3XxX, 



dN 
dY 

dN 

dZ Z r„ ,^ 9T 

dN - 2^^"^^ + " ]' 
— = - [3wxX + u -1]. 

Such a set of equations, together with the constraint 



(8) 



X + Y- 



(9) 



imposed by eq. ([2]), form the autonomous system upon 
which the remainder of this paper will be based. Two 
remarks are in order here. First, we define X and Y and 
not, say, x'^ and as variables to naturally allow for 
negative px and pc- While px < not only cannot be 
ruled out but in fact is mandatory in many /(i?) theories 
that aim to solve the coincidence problem and satisfy cur- 
rent observational constraints [23j, one could argue that 
Pc < is unphysical or at least undesirable. Nevertheless 
this class of models demands such flexibility, and to deny 
it a priori leads to an incomplete understanding of the dy- 
namics involved.'^ Any such physical considerations must 
be done a posteriori, for instance, by ruling out trajecto- 
ries for which pc < at some point. Second, due to the 
above constraint the system reduces to a 3— dimensional 



^ A previous work 1241 neglected this, and in so doing overlooked 
the existence of some fixed points and allowed for trajectories 
which would eventually violate the condition (|9jl. 



one, and therefore one of the equations (|8| may be elimi- 
nated. By doing so the reduced system still has to satisfy 
a constraint which, choosing to eliminate z, is written as 



X + Y- 



< 1. 



(10) 



Our notation for the fixed points will be in the form 
(ATfp, Ifp, Zfp, Ufp). From the last two equations of sys- 
tem Q it is straightforward to see that, for any fixed 
point, Z[p implies Ufp ~ Xfp = . This in turn 
forces either Yfp = or A^ = . The former option cor- 
responds to the baryonic point (0, 0, 1, 0) . The latter, to 
a regular CDM universe which is known to converge to 
(0, Ifp, Zip, 0) . All other fixed points will have Zfp = . 
On the other hand, when U[p ^ it can be shown that 
the only possible fixed point is (0, 0, 0, 1) , unless A^;, Ac 
and Wx are chosen to match a very specific relation given 

by 



(l~3wx)(l + 3Ac) =3A,. 



(11) 



Thus, the fixed points of interest, responsible for the 
MDE and present acceleration, can be written as 
{Xip, Yfp, 0, 0) . As mentioned before, we shall follow the 
approach of [8 and seek two such points: a saddle point 
responsible for the MDE, and an attractor governing the 
dark energy dominated epoch. 

These two important points will be dubbed 
A = (Xa,Fa,0,0) and B = (ATs, y^, 0, 0) . From 
the system ([S]) we find 



Xa 



Wx + Xx - Xc+ ^ {Wx + Aa; - Ac)2 ^WxXc 



(12) 



Xn = 



-I- Aj; - Ac - ^ (wx + Xx - Ac)2 + AWxXc 



2Wr, 



(13) 

with Ya,b = 1 — Xa,b ■ Note that Xb > Xa, since 
Wx < 0. In section IV we will show that points A and 
B are candidates for the MDE and present accelerated 
epoch, respectively. For the moment, it is useful to in- 
troduce the inequalities 



-Wx ^ \wx\ > (^V>^ + , (14) 

'Wx = \wx\ < - \/K:^ ■ (15) 



Note that the second condition on Wx is very strong since 
an eventual acceleration requires either A^; or Ac to be 
order unity for Wx to be negative enough, which is at odds 
with what we expect from previous results |24j . Both 
points exist whenever: (i) A^; > 0, Ac > and either 
condition (14 1 or (151 is met; or (ii) at least one of A^; 
and Ac is negative. However, depending on the sign of 
these parameters we can have Xa,b negative or larger 
than 1. If Xx and Ac are positive and condition (14 1 



is met, both points are placed between and 1. On the 



other hand if the inequality ( 15 1 is the one satisfied, both 
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Point 


X 


Y 


z 


u 


Existence 


Stability 


Acceler. 




Xa<0 


1-Xa 








Ac < 
[or] 

Ac > and A^ > Ac and < (\/A^ — ^A^) 


saddle under 
condition ( 16 1 
unstable otherwise 


no 


A2 


0<Xa<1 


1-Xa 








Aa; < and Ac > 
A:r,Ac>0 and \wx\ > {VX^ + VX^) 


saddle point 


WxXa < - | 




Xa > 1 


1-Xa 








< \x < \c and \wx\ < (\/A^ - \^)^ 


saddle point 


WxXa < - | 


Hi 


Xb < 


1 — Xb 








Ac > and \x > Ac and \wx\ < (\/Xe — \/X;)^ 


saddle point 


no 


B2 


< Xfl < 1 


1-Xb 








A:r > and Ac < 
A^,Ac>0 and \wx\ > {^/K + VX) 


attractor 


WxXb < -| 




Xb > 1 


1-Xb 








\x < 

< Aj; < Ac and \wx\ < [V-^x - vAcj 


attractor 


WxXb < -| 


C 








1 





Wx, Xx, Xc 


saddle point 


no 


D 











1 


Wx, Xx, \c 


unstable or saddle 


no 


E 





Ye 


ze 





Ac = 


saddle for > \x 
attractor otherwise 


no 


F 


—00 


+00 


+00 


+00 


undetermined, but see figure 3 


N/A 


N/A 



TABLE I: The properties of critical points for the Dark Interactions Model. It is implicitly assumed that Wx < and that 
condition (11 1 is not met. The last column show the conditions for late-time acceleration of the universe. Figure [l] illustrates 
this table. 



Xa and Xb will be negative (for > Xc) or larger than 
one (for A^; < Ac). When A^: < , we have Xb > 1 and 
sign(X^) = sign(Ac). Lastly, in the case Ac < , we find 
that Xa < and Xb < 1 if and only if A^; > 0. 

Another possibility is the existence of a catastrophical 
"point" at infinity (X —00, Y —^ +00), to which 
the system can collapse in some cases (the reciprocal 
point with positive X and negative Y can be shown not 
to exist). To illustrate such a case let's consider, for 
the sake of argument, Xx < , Ac>0, Y > and 
z ~ u ~ . Then, the term — 3Ac Y is negative and 
we see from (|8| that dX/dN will be negative for small 
enough positive X, and will stay so after X becomes neg- 
ative. At the same time, dY/dN will be positive for large 
enough values of Ac. A catastrophic decay thus ensues 
and X —00 in finite N."^ 

Finally, as pointed out before, the system has three 
other points of interest. A baryonic point C, a radia- 
tion dominated point D, and a mixed matter dominated 
point E, defined respectively as (0,0,1,0), (0,0,0,1) 
and (0, Y, z, 0) . We find that point C is always a sad- 



In this case none of the energy densities are actually growing to 
infinity, but rather it is H that is going to zero. 



die point, point D may be either unstable or a saddle 
point depending on the values of Xx , Ac and Wx , and that 
point ii' is a saddle point unless \wx\ < Xx, in which case 
it is an attractor. 



A. Stability of the fixed points 

Table |T] summarizes the conditions for existence and 
stability of all 6 fixed points. For the sake of clarity, we 
subdivide each point A and B into three, according to 
whether Xfp < 0, < Xfp < 1 or Xfp > 1. 

Regardless of the positivity of A^; and Ac, points A2, 
A3 and Bi are always saddle points, whereas points B2 
and B3 are always attractors. Point Ai is a little trickier. 
It will be a saddle for either 



|w^|(l-f 3Ac) > A^-Ac- 



[Ax — Ac 



1 



(1 + 3Ac) 



< \Wx\ < 



9Ac 



if Ac < ; 
if Ac > . 



(16) 



In all other cases, it will be unstable. In any case, whereas 
points A2, A3 and Bi have 2 negative (and 1 positive) 
eigenvalues, point Ai may have at most only one negative 
eigenvalue. 
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FIG. 1: Phase space trajectories together with the fixed points 
A, B, C and D. From top to bottom we have Aa: = Ac = , 
Ai = Ac = 0.04 and A^ = —Ac — —0.15, while Wx ~ —1 
for all three plots. Note that some trajectories cross either 
the X = or Y = surfaces. The gray volume represents 
the region where the constraint ( |10[ | is violated. The black 
dashed fine corresponds to a trajectory that passes through 
the present (observational) energy densities, which in turn is 
depicted by the red star. 



The baryonic point C is a saddle point, but since we 
do not expect a baryonic dominated era, for reasonable 
initial conditions trajec tories should not pass too close to 
it (although see section VII I . As for the radiation domi- 
nated point D, even though it may technically be a saddle 
point (depending on the values of A^;, Ac and Wx), in all 
physically relevant scenarios it will be unstable. 



B. Summary of fixed points 

Throughout this subsection we shall assume, for 
simplicity, that neither peculiar conditions (111 nor 
Xc — hold. We thus have a total of 4 fixed points plus 
a possible "catastrophical abyss" dubbed F. Point D 
is the initial condition for our system. Proximity to 



point C is better avoided in order to cope with obser- 
vations. Points A and B are the main focus of interest: 
the former, a candidate for a saddle MDE and the latter, 
a candidate for the present accelerated epoch. Point A 
and B will usually be a saddle and an attractor, respec- 
tively, but there are exceptions. 

In this model, contrary to the particular case |24j . we 



have 



Ya 



Xb- 



In fact, we have Xa — Yb 



Aa; — Ac 



instead of zero. This is a desirable result, for it allows 
the accelerated fixed point to have a dark energy density 
closer to the present value. Put in another way, when 
one adjusts the MDE fixed point found in [53] to cope 
with observations, the accelerated point acquires a ratio 
^ = Pa;/Pc ^ 45 (see sections |v] and VII I, while present 
observations state that Rq k, 3.4. We will come back to 
this issue on section |VII| where we investigate possible 
improvements when the two coupling constants differ. 

A sample of possible trajectories are depicted in fig- 
ure [l] for three different cases, all of which have Wx set 
to -1: A^ = Ac = (ACDM), A^ = Ac = 0.04 and 
Aa; — —Ac — 0.15. The gray volume denotes the forbid- 
den region, i.e., the one for which the constraint (10 1 is 
violated. The black dashed curve represents a trajectory 
that passes through the present (observational) energy 
densities, which in turn is depicted by a red star. Note 
that, as pointed out before, in all cases it is possible for 
trajectories to cross either the X — Q or Y — Q sur- 
faces. 



Had we chosen 



as variables, we would 



' Px/e 

see trajectories hitting such surfaces but artificially not 
being able to cross them. This is specially important for 
models with negative values of A^; and/or Ac, for which 
the fixed points A and B may have either Xf-p < or 
Ytp < 0. 



IV. POSSIBILITY OF TWO SCALING EPOCHS 

We now analyse the feasibility of a duo-scaling cos- 
mology, in which the MDE is assured by the existence 
of a suitable saddle point and the present acceleration 
by the existence of an appropriate final attractor. From 
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table [T] it is clear that we need point A to be such a sad- 
dle, with Xa ^ 0.1, and point B to be the attractor, 
with Xg > 0.7 . In other words, points and Bi are 
excluded. Point B^, in turn, is not very interesting for it 
would mean that: (i) the universe is now in a transient 
phase between the two scaling epochs, which inevitably 
requires a more precise tuning of the parameters; and 
(ii) the dark matter energy density would become nega- 
tive in the future. The former is a crucial feature, as it 
automatically prevents solution to the CP in the way we 
defined it. We shall thus consider only the B2 flavour of 
point B and consequently require Xx > 0. In a nutshell, 
we limited ourselves to only two possibilities: 

• Point Ai followed by B2 (requiring A,, < 0); 

• Point A2 followed by B2 (requiring Ac > 0). 

To make point Ai a viable saddle we need Ac to be 
negative, since conditions (14 1 and (15 1 are mutually ex- 
clusive. Also, from eq. (16 1 and the fact that A^; > , 
we find that Ac must be larger than —1/3. But even 
— 1/3 < Ac < might not be enough, since point Ai has 
at most one negative eigenvalue and, as such, it is unclear 
if it would have a large basin of attraction as required to 
address either the fine-tuning of initial conditions (SIC) 
or the CP. In other words, the corresponding MDE might 
be short and it may lead to smaller values of A. 

Points A2 and B2 are easier to reconcile. In fact, a 
sufficient condition for the system to have the two scaling 
regimes we need in other to attack the coincidence and 
initial condition problems is A^; , Ac > and Wx satisfying 
the constraint (14 1. However, as will be shown in the 



next section, observational constraints put limits on the 
coupling strength which are tighter in this case than in 
the one with Ac < 0. 

Following these guidelines, we develop a quantitative 
analysis in section |VT] for some of the observationally al- 
lowed values of the parameter space. We also investi- 
gate the possible improvements on the range of initial 
conditions that eventually evolves to today's observable 
universe (i.e., the DEICP). 



V. OBSERVATIONAL CONSTRAINTS 

In order to probe quantitatively the permissible range 
of values of the two coupling constants, we submit the 
model to two different observational tests: type la su- 
pernovae and the so-called CMB shift parameter. The 
former relies on the well accepted hypothesis that this 
kind of astronomical object is a standardizable candle, 
and consists of comparing their distance moduli at dif- 
ferent redshifts to the ones calculated from the model. 
The latter can be considered a good first approximation 
of the full CMB analysis (specially when used in con- 
junction with the acoustic peak scale [26]), which would 
require setting up the first order perturbation equations 
for the model and employing one of the established CMB 
codes, such as CAMB or CMBEASY. 



We made use of the combined Essence, Hubble, SNLS 
and nearby supernovae catalog as compiled by [lOj^, for 
a total of 192 supernovae, and of the WMAP 3-year re- 
sult for the shift parameter [TT] (to wit, 1.70 ± 0.03). 
We employed a grid-based method to compute the for 
different values of A^;, Ac, Wx and f2cO- The present val- 
ues of the baryonic energy density fl^Q and the radiation 
energy density ilro were held fixed respectively at 0.042 
and 4.2 x 10~^, which are the best fit values of the com- 
bined WMAP3 and SDSS observations |27]. The value 
of h was marginalized analytically. For f2cO we assumed 
a gaussian prior equivalent to WMAP3-f SDSS observa- 
tions, i.e., with a mean at 0.22 and a standard deviation 
of 0.03. We also considered a top-hat prior in the range 
0.14 < fico < 0.30, but results were very similar. 

Figure |2] depicts the 1 and 2a confidence levels for the 
gaussian prior and contains much information about the 
model. It shows that, in absolute terms, negative values 
of Xx can be many times larger than positive ones. It is 
also clear that neither test disfavours phantom values of 
Wx, and that the allowed region in the A^^ > 0, Ac > 
quadrant is quite small. A two dimensional "cut" at 
Wx — ~l is shown in figure [3] and gives a good idea of the 
range of values A^; and Ac can take as well as the overall 
alignment of the contours with the line A^; — —kXc where 
fc ~ 2. Technically it is not a cut from figure [2] but rather 
the contours for a two parameter model with Wx set at 
— 1 from the beginning. In practice the only difference is 
that the contours in figure [3] are smaller than the yellow 
ones in figure [2j Figure [4]shows different two dimensional 
contours derived from the higher dimensional likelihood, 
assuming flat priors on all 3 model parameters, as well 
as the one dimensional likelihood for each one. The use 
of flat priors allows the interpretation of all these six 
plots as projections of the higher dimensional ones, thus 
helping visualize figure |2] Finally, figure [5] merges the 
three dimensional plot with the regions of the parameter 
space that exhibit the duo-scaling regime. In the inte- 
rior of the green transparent region we have a, Ai — B2 
duo-scaling, whereas in the interior of the gray opaque 
checkered volume we have a, A2 — B2 one. In addition to 
the existence and stability of the fixed points A2 and B2 , 
the gray borders are drawn in such a way to guarantee 
that Xa2 < 0.1. The idea is that dark energy needs to be 
sub-dominant during the MDE in order not to interfere 
too much with the formation of structures. In any case, 
as can be seen from the plot, a fixed point Xa2 > 0.1 is 
excluded at over 2a. 

The major focus of previous works which adopted sim- 
ilar coupling models has been in the positive A^;, Ac quad- 
rant. It is now clear that such a region is a severely lim- 



^ This catalog has the weU-known feature of agreeing with ACDM 
within 1(T. Here we are not advocating its use over other SNIa 
sets, but the reader should bear this in mind when analysing our 
confidence levels contours; to wit, the inclusion of interactions 
do not change this picture. 
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FIG. 2: Combined results from supernovae and shift param- 
eter tests. The volumes represent 1 and 2a confidence lev- 
els, marginalized over fico with a gaussian prior based on 
WMAP3-I-SDSS results. The vertical contours are drawn at 
\x = ^ (white) and Ac = (black). The yellow horizontal 
cut is made at = — 1 (see figure [sj , and further dashed 
cuts are made to aid the eye at = —1.5, = —2.0 and 
~ —2.5. 

ited piece of the allowed parameter space. From figure [3] 
we see that if we limit ourselves to, say, A^; = Ac > 
as done in 25], we must have Ac < 0.014 (0.020) at 
2(7 (3(t) confidence levels. If we consider negative Ac, on 
the other hand, the coupling constants can be as large 
as 0.30 (0.42) and -0.20 (-0.35) at 2cr (3cr) for A:^ and 
Ac, respectively. That is over an order of magnitude in- 
crease on absolute values. We will exploit such freedom 
in what follows and investigate the feasibility of models 
with negative values of Ac. 



VI. FINE-TUNING AND THE INITIAL 
CONDITIONS PROBLEMS 

As pointed out before, there are two ways in which any 
model that aims to address the DEICP with small SIC 
may offer improvements. One is to increase the value of 
the its initial DE energy density at some given time in the 
past, or more precisely, to increase the ratio C of energy 
densities between dark energy and radiation at the end 
on inflation. The other is to increase the range A which 
will evolve to the observed values of ^xQ- The former 
may help conciliate the dark energy component, under 
a field interpretation, with our present understanding of 
the Standard Model of particle physics or of its exten- 
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FIG. 3: 1, 2 and Sa confidence level contours for lu^, fixed 
at —1. The brown dashed line shows the particular case of 
= Ac , as considered in [T7], and the yellow dot stands for 
the ACDM case. The gray area on top is the region for which 
the system enters the "catastrophical abyss" characterized by 
the fixed point F, which in practice means that H{z*) — 

for some Z <C ^recombination • 



sions. The latter is a direct measure of the basin of at- 
traction of the MDE saddle point and, therefore, of the 
reduction of fine-tuning required. 

A model which undergoes a duo-scaling regime is ex- 
pected to be less sensitive to initial conditions due to the 
presence of the MDE saddle point A. The amount of 
loosening is nonetheless unclear from the previous argu- 
ments. One needs therefore to assess the range of allowed 
initial conditions for, say, the energy density of dark en- 
ergy for different values of the parameter space allowed 
by observations and compare it with the range of initial 
conditions in the ACDM model. This is done in figure |6] 
where for = — 1 we plot the values of A/Aacdm 
and C by computing the range of initial values of px that 
evolve to present day observations {Xq — 0.74 ± 0.03). 
As noted in section |l] we take our initial values at the 
end of inflation, more precisely at z = 10^^. For refer- 
ence, Aacdm — 0.08 (i.e., initial conditions for may 
vary ~ 8%, which is of course what we expect since in 
ACDM px is constant and that is its current uncertainty) 
and Cacdm — 9 x 10"^"^. We see that when interactions 
are at work, improvements of a factor between 1 and 30 
for A and of dozens of orders of magnitude for C (as 
long as Ac 0) are possible. For small but positive val- 
ues of Ac (starting around 10""*), ( drops abruptly to the 
ACDM value. We remark that the negative Ac region has 
negative pxi, and thus in this section ^ looses its impor- 
tance since the DEICP is no longer applicable (remember 
that we are identifying the DEICP with equipartition ar- 
guments) . Note that the valley close to Ac = , and 
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the thin peak at its base (A^ ~ —0.18) in the plot of A 
were artificially widened in figure |6] for clarity purposes. 
Both are actually much thinner than depicted, and only 
exist for |Ac| < 10^^". As a matter of fact, this peak 
at the bottom is spurious: p™^'', the denominator of A, 
is getting very close to zero there. As predicted in sec- 
tion |IV| the region of positive A^; > is less sensitive 
to initial conditions, which is due to the presence of the 
duo-scaling regime. 

Two important remarks are in order here. First, as 
far as the DEICP is concerned, it is much more sensi- 
tive to Ac than A^,. In fact, the highest possible values of 
C are obtained for very small and positive Ac. Stronger 
couplings with DM only make the DEICP worse, where 
as no values of A^; allowed by the observational tests ap- 
plied make a difference when Ac ^ 1. Second, in the 
region of positive A^;, which is the interesting one since 
in allows the existence of both scaling regimes, the SIC is 
considerably smaller than in ACDM, which is an advan- 
tage for any phenomenological approach. In particular. 



the region around A^; = 0.2 and A^; — —0.1, which show 
improvements in A by factors of 10 — 30, turn out to 
be the most promising one to solve the CP, as will be 
discussed in the next section. 



VII. A TWO- VARIABLE DUO-SCALING TOY 
MODEL 

Inspired by the results on the preceding sections, we 
look for a toy model which at the same time alleviates 
the DEICP, do not has a large SIC, undergoes a duo- 
scaling regime and satisfies the observational tests carried 
out in this work. A good starting point is to select a 
model which has an analytic solution. As mentioned in 
section |TT] the best candidate is the one in which 

Wx = \x ■ (17) 
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FIG. 5: Intersection of the 1 and 2a contours with the regions 
of the parameter space that aUow for a duo-scahng regime. 
Inside the green transparent region we have a. A\ — B2 type 
of scaling, whereas inside the gray opaque checkered region 
the scaling is of the type A2 — -B2. The gray borders are 
drawn in such a way to guarantee that Xa.^ < 0.1, so as not 
to compromise the formation of structures during the MDE. 
The 2D contours are the same as in figure [2] 



In this case, we can write the ratio R between and pc 
as 



Logiol^] 




FIG. 6: The DEICP variable C (top) and SIC variable A 
(bottom), plotted as a function of the coupling constants Xx 
and Ac. Here we set Wx ~ —1 ■ The improvements in range 



between 55 and 75 orders of magnitude (Cacdm 



10" 



0, 



while A is between 1 and 30 times the ACDM one. The abrupt 
decline of ( close to = is expected (see section VII I. The 
width of the valley at Ac = in the plot of A is grossly 
exaggerated here for clarity purposes. In fact, the drop only 
starts at |Ac| < 10"^". For Ac < 0, is negative, and thus 
looses its relevance since the DEICP is no longer applicable. 



R 



jS + L)Ro - 2A, + (1 + 0)3-^ [{S - L)Ro + 2X^] 

S - L + 2RaXx + {l + z)3S [{S + L) - 2i?o A^^] ' 

(18) 



where S = ^/Y- 

Ro = Pxo/PcO - 



Ac(-2 - 4Aa; + Ac) , L = 1 - Ac and 
3.4 . To better understand the in- 
fluence of the two couphng constants it is fruitful to 
write down the above expression in hmits z — > cxd and 
z —^ —I. Since we will only be considering coupling con- 
stants which are reasonably smaller than unity, we will 
neglect all but the two lowest order terms in Aa; and Ac 



in (18). We thus get 



lim i? = Ac(l + Ac-i?oA:r) + 0(A3), (19) 

; — *oo 

Ro + K - Ac(l + i?o) 



lim R 

z^-l 



Aa; 



+ 0(A) . (20) 



We see that for Ac = 0, R goes asymptotically to zero 
at high redshifts. If Ac in nonzero we get (unless very 
specific values of the couplings are chosen to cancel the 
term proportional to 1 -I- z in the numerator) a constant, 
nonzero value for R as z goes to infinity. On the other 
extremity, as z — > — 1, we see that setting A^ = means 



that the ratio of energy densities will grow indefinitely in 
the future. This gives a physical intuition on the roles of 
both Xx and Ac when they assume nonzero values. The 
former is needed to ensure a late-time scaling, i.e., to 
keep R small, while the latter ensures a nonzero value of 
R at early times. Put in another way, large values of A^; 
help solve the CP, while large values of Ac alleviate the 
DEICP. Referring back to section |III| this could also be 
seen from the fact that = when Ac = and that 
= 1 when Aa; = 0. 
We now need to select particular values for Aa: and Ac 
(remember that we are assuming Wx = —1 — Aa:) that 
passes the proposed observational tests and gives a fi- 
nal value R^i of R which is not far from the present 
one. From figures [2] and [3] we see that the lowest asymp- 
totic value of R allowed in the positive coupling con- 
stants quadrant is ~ 32 (i?_i ~ 45 if Aa: = Ac), 
while for negative Ac we can get as low as ~ 5. 

In any case, the coincidence problem is at best only 
alleviated, not solved, as we would still be some- 
what far from these values. Based on these consid- 
erations, we select {Aa; = 0.18, Ac = —0.08, Wx = —1.18} 
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FIG. 7: Evolution of all four energy densities as a function 
of redshift. From left to right, the dominant component is, 
in order: radiation, dark matter and dark energy. Baryons 
are always sub-dominant. The dashed lines correspond to the 
ACDM model and the solid lines to the case \x ~ 0.18, 
Ac = — 0.08 and w; a; = —1.18. Note that matter-radiation 
equality is pushed back to 2 ~ 20000, and that DE has neg- 
ative energy density for 2 > 2. 
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FIG. 9: Same as figure [7] for the case where = 0.23, 
Ac = —0.14, Wx ~ —0.85. Again, matter-radiation equal- 
ity is pushed back to 2 ~ 20000. Note that for a not-so-short 
redshift interval, baryons become the (marginally) dominant 
component. In this scenario, we have already reached the ac- 
celerated attractor (R-i within 4% of Rq) and thus, by our 
definition, solved the CP. 



10'' 10'' 10^ 10- 



0.75 



G 0.5 



0.25 



0.75 




0.25 



FIG. 8: Same as figure[7]for the case where Ai = Ac = 0.014, 
Wx = —1.014. Here the evolution resembles more closely that 
of ACDM, and in particular the DE energy density is always 
positive. 



and {Xx — Xc — 0.014, Wx = —1.014} as examples of our 
toy model. In the former, early DE has negative energy 
density, and thus it is not a candidate for solving the 
DEICP. The latter has positive dark energy throughout 
the whole history, but exhibits a worse CP. 

In figure [7] we plot the diff'erent energy densities as a 
function of redshift for the first chosen set of parameters 
and compare them with ACDM. The present values of the 
difi'erent are set by the WMAP3-fSDSS best fit. Note 
that px is negative throughout the entire matter domi- 
nated era, which helps to increase its duration. Another 
important feature of this example is that the radiation- 
matter equipartition occurs much earlier, at z ~ 20000, 
and therefore it is possible that different observational 
tests such as a full CMB analysis would rule out such 
an anticipation of the MDE. In fact, fitting the observed 



matter power spectrum would be a challenge since its 
peak is roughly estimated by fcoq, the wave number of per- 
turbations which enter the horizon at matter-radiation 
equality, and that would be shifted to smaller scales by 
a factor of ~ 6.1 (in units of ft. x distance"""^). This 
example has small SIC, as can be seen from figure |6] 
Figure [8] is the same for the second example, for which 
Xx — Xc — 0.014. In this case, the DEICP does apply 
and we get C — 3 x 10~^^ which should be compared 
to Cacdm — 9 X 10~^°^. Note that in this case we are 
farther from the final accelerated attractor than in the 
previous one. In fact, R-i = 6.1 in the first example and 
= 70.4 in the second one. 
We conclude this section by noting that it seems to 
be possible to get R-i closer to Rq for regions of the 
parameter space that do not allow analytic solutions. 
Inspired by eq. (20 1, we seek large values of Xx, and 
from figure [2] we see that larger A^; are allowed when 
Wx > —1- Probing that region, we found a good candi- 
date for {A^ = 0.23, Ac = -0.14, Wx = -0.85}. Figure [9] 
depicts the evolution of the model for these parameter 
values, which give R-i = 3.48. This is within 4% of 
the present (WMAP3-I-SDSS) value, and one could ar- 
gue that this solves the CP. Nevertheless, once again the 
early equality might jeopardize the model (here fceq is 
shifted by a factor ~ 6.6). Another aspect we must point 
out is that this region of the parameter space is a rather 
sensitive one. This is refiected in the wiggles and dis- 
continuities on the tip of the contours of figure |2] What 
this means is that small changes in the parameter values 
are "amplified", i.e., may result in substantially different 
evolutions. For instance, for A^: — 0.235 dark matter 
has negative energy density for z > 8; for A^; = 0.20 
the matter-radiation equality is pushed back to z > 10^. 
This can be problematic if one tries to extend our anal- 
ysis for varying Wx and/or varying Aa;, Ac scenarios. 
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VIII. CONCLUSIONS 

In this paper we tried to clarify some of the different 
aspects of cosmological fine-tuning and their relation to 
dark interactions. With this goal in mind, we proposed 
the use of three variables to quantify each of most com- 
mon issues of cosmological models: the CP, the DEICP 
and the SIC. By restricting our parameters through the 
use of both supernovae and the CMB shift parameter and 
through the requirement of a duo-scaling cosmology, we 
greatly limited our parameter space. 

Applying the proposed variables to our Dark Interac- 
tions Model, which has some distinct features such as 
negative DE energy densities, we found that each cou- 
pling constant is related to a different regime. The CP 
can only be solved for large values of A^^, while larger 
values of Ac guarantee higher amounts of DE in the early 
universe, and thus relate to the DEICP. We also found 
that nonzero values of Ac also give rise to smaller SIC for 
any value of Aj.. We thus investigated two examples of 
a toy model class (characterized by eq. (17l) for which 



there exists analytic solutions so as to gain some intu- 
ition on how to address these different problems. For 
both of them the CP and DEICP were only at best alle- 
viated, since the lower value of R^i obtained was 6.1 and 
the higher value of C was ~ 10^^^. Even if not a proper 
solution, the latter represent a 74 order-of-magnitude im- 
provement over the ACDM scenario. It is possible that a 
solution to the DEICP requires a third scaling solution: a 
scaling between dark energy and radiation. Such regime 
would only be achieved if we shed light on the dark in- 
teractions, including radiation in our coupling scheme. 



It is also noteworthy the fact that, based on the exam- 
ples shown, solutions to the CP (in our class of models) 
seem to require DE with negative energy density in the 
past. This may be a hint that modified gravity theories, 
which can more easily accommodate negative (effective) 
energy densities, might have an upper hand at explain- 
ing this issue. In fact negative-energy fields spreading 
over large spatial regions might yield a number of exotic 
phenomena (see |28j and references therein). 

Finally, we concluded that ~ Rq only occurs in 
a rather turbulent region of the parameter space, for 
which Wx is not very negative and A^, is large enough. 
Since in this case the dark matter-radiation equality dif- 
fers substantially from that of ACDM, it remains to be 
seen if further observational tests (specially those involv- 
ing perturbations) will rule it out. In particular, coping 
with the shift in the peak of the matter power spectrum 
will be a challenge for this model. Should the perturba- 
tion equations turn out to be considerably different than 
those of the concordance model and the matter power 
spectrum fitted, then as sensitive as the model might be 
to its parameters in this region, compared to ACDM or 
quintessence no real fine-tuning will be necessary in order 
to solve the coincidence problem. 
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